Coupled Superconducting Phase and Ferromagnetic Order Parameter Dynamics 
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Via a direct coupling between the magnetic order parameter and the singlet Josephson supercur- 
rent, we detect spin- wave resonances, and their dispersion, in ferromagnetic Josephson junctions in 
which the usual insulating or metallic barrier is replaced with a weak ferromagnet. The coupling 
arises within the Fraunhofer interferential description of the Josephson effect, because the magnetic 
layer acts as a time dependent phase plate. A spin- wave resonance at a frequency u s implies a 
dissipation that is reflected as a depression in the current- volt age curve of the Josephson junction 
when %ujs — 2eV. We have thereby performed a resonance experiment on only 10 7 Ni atoms. 



The coupled dynamics of the electromagnetic field and 
a Josephson junction has a number of manifestations 
and is very well understood [l], 0, 0, 0] • When the 
usual insulating or metallic barrier is replaced with a 
weak ferromagnet there is a coupling to another field, 
namely the spontaneous magnetisation of the ferromag- 
net. Spin- waves are elementary spin excitations which 
can be viewed as both spatial and time dependent vari- 
ations of the magnetisation. In a ferromagnet the lowest 
energy excitation, the Ferromagnetic Resonance (FMR), 
corresponds to the uniform precession of the magnetisa- 
tion around an externally applied magnetic field at the 
frequency uo s . This mode can be resonantly excited by 
an alternative (ac) magnetic field that couples directly to 
the magnetisation, as described by the Landau-Lifshitz 
equations [1] . The Josephson phase difference (j) between 
the two superconductors has its own dynamics. A bias 
voltage Vo causes qb to become time-dependent so that 
</> = </>o + ujt, where uoj = (2e/h)Vo and qbo is arbi- 
trary. Corresponding to the ac Josephson effect [l[, for 
our junctions, to a good approximation, the resulting ac 
Josephson current density is J s = J c sin((/>o +c<;jt), where 
J c is the critical current density. 

In analogy with the A-phase @ of 3 He, coupled mag- 
netic and phase oscillations should exist in ferromagnetic 
superconductors with triplet pairing, but have never been 
observed. We show here that a similar coupling for sin- 
glet superconductors can be realised in a Josephson junc- 
tion with a ferromagnetic barrier. The dynamical cou- 
pling stems from the spatial interference of the Aharonov- 
Bohm phase caused by M(£), resulting in the spatial de- 
pendence of 0(r,t). The ac Josephson current produces 
an oscillating magnetic field H(t) and when the Joseph- 
son frequency matches the spin wave frequency, uoj ~ uo Sl 
this resonantly excites M(t). Due to the nonlinearity of 
the Josephson effect, there is a rectification of current 
across the junction, resulting in a dip in the average dc 



component of J s at voltage V s = (h/2e) uo s . The principal 
result reported here is the observation of these coupled 
dynamics. 

Magnetised Josephson junctions Q require weak ferro- 
magnetic materials and nanosized junction area to keep 
the overall magnetic flux in the junction smaller than the 
flux quantum $0. An electron microscope image of a typ- 
ical ferromagnetic junction used in this study is shown in 
Fig. 1(a), while Fig. 1(b) is a schematic representation of 
the different layers. The superconducting electrodes com- 
prise 50 nm of Nb (T c = 7.6 K), while the barrier is 20 nm 
of Pd . 9 Ni .i (Tcurie = 150 K). The current-voltage (IV) 
characteristics are measured using current bias and are 
reported, as function of the applied in-plane field, in the 
right insert of Fig. 2. The IV characteristics are not hys- 
teretic, and overall they correspond closely to those ex- 
pected for a junction with a conductive barrier [1, 0]. 
The junction normal resistance is R n « 0.8 ^, and the 
Josephson coupling is I c R n ~5 /iV, as expected for ferro- 
magnetic junctions of this thickness [7], yielding the criti- 
cal supercurrent of I c « 7 fiA. The hostile nature of even a 
weak ferromagnetic environment for singlet Cooper pairs 
is illustrated by a similar junction with 70 nm of non- 
magnetic Pd which, despite the almost four times larger 
thickness, has a larger critical current I c « 44 juA. 

For a square junction of side L, the total supercurrent 
is given by the integral [To[ 



fL/2 rL/2 

I s = J c dx dy sin <j>(x,y,t) (1) 

J-L/2 J-L/2 



with 



2e 

(r, t) = 0o + u jt - — j A • dr, 



(2) 



where the last term is the Aharonov-Bohm phase [111 ], 
involving the vector potential A. We use a gauge where 
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FIG. 1: (color online), (a) SEM photo of the ferromagnetic 
Josephson junction used in the experiment. Junction area is 
500 nm x 500 nm. (b) Schematic cross-section of the Joseph- 
son junction. The ac Josephson current I(t) flows through the 
junction creating an rf magnetic field H(t), causing the mag- 
netisation precession M(£), which in turn resonantly couples 
with the Josephson phase at frequency coj. Layers are respec- 
tively from the bottom: Nb (50 nm), Pdo.gNio.i (20 nm) and 
Nb (50 nm). (c) The equivalent RSJ model and the sketch of 
the effect of the FMR on the current- volt age characteristics. 
The ferromagnetic layer is modeled as a series LCRo oscilla- 
tor, in parallel with the Josephson junction. Resistance Ro is 
proportional to the imaginary part of the susceptibility x" • 



A = A(r, t) z, the direction z being perpendicular to the 
junction surface [see Fig. 1]. Therefore 0(r, t) = </>o + 
kx + uojt + </> m , where (j) m = (4ae/K)A mz reflects time 
dependent fields and k = (4ed/Tt)/jLoH + (Aea/fijfioMoy. 
Here Mq v is the y component of the static magnetisation 
M , the applied field H is in the y direction, 2a and 
2d = 2 (a + A) are the actual and magnetic thickness of 
barrier and A the London penetration depth. Equations 
(PQ) and ([2]) are used to describe both the statics and the 
dynamics of our junctions. 

In the absence of a bias (Vb = 0), we are dealing 
with static fields, and the Equations (pQ) and ([2]) lead 

rL/2 



, dx sin kx [11]. 



to the Fraunhofer pattern I s = J C L J_ L j 2 
The magnetisation Mq of the barrier has the same effect 
as inserting a wedge shaped phase plate in front of the 
slit, it displaces the diffraction pattern. Experimentally, 
the diffraction pattern is shifted to the right for increas- 
ing (positive Mq v ), and the left for decreasing (negative 
Mq v ), fields. This illustrates the linear nature of the 
coupling to M. In Fig. 2, the dotted curves are a fit us- 
ing Eqs. (pQ) and ([2]), along with the magnetisation data 
measured on a trilayer with the same cross section as the 
junction [see the left insert of Fig. 2). The periodicity 
and the asymptotic behaviour of the measured diffrac- 
tion pattern attest to the high quality of our junctions. 
They confirm the close-to-uniform current distribution 
and single-harmonic current-phase relation, while the re- 
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FIG. 2: (color online) Dependence of the Josephson critical 
current on the external magnetic field in plane. The solid 
curve represents the normalized experimental data taken at 
35 mK and the dotted curve is the Fraunhofer pattern ex- 
pected for our junction parameters including the magnetisa- 
tion. Left insert: The magnetic hysteresis taken at 10 K on 
a millimetric trilayer with the same cross section as the junc- 
tion, with field applied in plane (dotted curve) and perpendic- 
ular to the plane (solid curve). Right insert: Current- volt age 
curves for different field values, measured at 35 mK. 



production of the shift with the two sweep directions, 
using experimental magnetostatic data, confirms the va- 
lidity of our description. 

The dynamical coupling reflects a similar phase con- 
tribution (/?(r, t) due to M(r, £), but which now has both 
a temporal and a spatial dependence, the equivalent of 
a phase plate in the optical analog with a similarly de- 
pendent refractive index n(r, t). The dynamics of the 
magnetisation (timescale of 1 ns) is much slower than 
the diffusion time through the ferromagnetic layer (0.5 
ps). The Josephson coupling is thus adiabatic with re- 
spect to the magnetisation dynamics. This assumption 
is implicit in Eqs. (pQ) and ([2]). The signal is seen for 
V > I C R, implying Eq. (pQ) can be linearized. The dc 
magnetic signal then corresponds to [Io[ 



4ae f L / 2 f L / 2 

I m = — / dx dyJ c cos(kx^uj J t)(l) rn , (3) 

n J-L/2 J-L/2 



where the bar denotes a time average. Substituting for 
mz and using J — V x H, following both 
a time and space integration by parts, the dc signal re- 
flecting the magnetic resonance is 



VoJ dt 



(4) 



with Mi(r,i) = / dt f Xi(t — t f )Hi(r,t f ), i = x,y,z, where 
Xi(t) is the dynamic susceptibility. This has an appealing 
interpretation in terms of magnetic losses. Here, as illus- 
trated by Fig. 1(b), H(r,t) is the magnetic field which 
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circulates inside the junction by virtue of the ac Joseph- 
son current. The junction lateral size L is smaller than 
both A and the skin depth for the frequencies involved. 
The displacement current is therefore negligible and all 
that is needed is to integrate Ampere's law in order to 
determine H(r, t). More details of these calculations are 
given elsewhere [121 ]. The current due to M(r, t) is 

I m = 2tt/ c (0)|^ [FMwj) + F vX '^j)\ , (5) 

where <£ r f = (2aL)B T { = (2aL){ioI c (0)/ L is the flux due 
to the radio frequency field and F x = (l/12)(7 c (£?o)// c ) 2 
and F y = (2/9 2 L )[l + sin((9 L /2) cos((9 L /2) + ((13/12) - 
(4/6>D)sin (# L /2)]; L = fcL, reflect the geometrical 
structure of the coupling. As the equilibrium magnetisa- 
tion is along the z axis, the magnetic resonance signal is 
contained in Xx( u j) anc ^ Xy&j), the Fourier transforms 
of the imaginary part of the susceptibility. Therefore, the 
total dc current within the Resistively Shunted Joseph- 
son junction (RSJ) model [1,0] is 



v foV) 

100 



V , I C 2 (B ) 



R(0) 



2V n 



R(loj) - I n 



(6) 



where R(0) and R(uoj) « R are the junction resistances 
for dc and frequency uoj. A simple physical argument 
can account for the three terms in Eq. (|6|). The average 
power dissipated in the junction is IVo and so the first 
term, Vq 2 /R, corresponds to the Ohmic loss at dc, while 
^I c 2 R(uj) is the similar loss at uj. The key third term 
represents a self-inductance L(M), stemming from the 
ferromagnet, in parallel with the junction and modeled 
as an LCRo oscillator [see Fig 1(c)], where R$ reflects 
the magnetic damping. At the magnetic resonance fre- 
quency, energy is absorbed by the ferromagnet, causing 
the oscillator to be lossy. This actually reduces the effec- 
tive junction resistance, leading to a dip in I(V). In this 
manner, the Josephson junction rectifies the self-induced 
magnetic resonance. 

This coupling to the magnetic system is evident in the 
measured dynamical resistance dV/ dl curves reported in 
Fig. 3. We measure the dynamical resistance rather than 
the IV characteristics to improve amplitude resolution. 
The mode labeled FMR (Ferromagnetic Resonance) is 
seen only for ferromagnetic junctions. There is good 
agreement between the experiment, solid curves, and the- 
ory, dotted curves. The magnetic resonance mode ob- 
served in our experiments reflects the properties of a thin 
film of the ferromagnet Pdo.9Nio.1- Magnetisation curves 
M(H), measured directly for a large area Nb/PdNi/Nb 
trilayer with the same cross section as the junction, are 
shown in the insert of Fig. 2. They indicate that M is 
perpendicular to the junction plane, a conclusion rein- 
forced by earlier anomalous Hall effect measurements on 
similar thin films [13]. The FMR mode, shown in Fig. 3, 
occurs at Vq = 23 /iV. This is unambiguously identified 
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FIG. 3: (color online). Dynamical resistance of the fer- 
romagnetic Josephson junction (solid curve, SFS, bottom 
and left axis) shows resonances compared to a similar non- 
ferromagnetic junction (solid curve, SNS, top and right axis). 
Dotted curve is a fit to theory [Eq. JBJ]. The mode at 23 jiV is 
the ferromagnetic resonance (FMR). Bottom insert: Conven- 
tional cavity ferromagnetic resonance on a macroscopic tri- 
layer. Top insert: Comparison between the field dependence 
of the FMR in the Josephson junction (solid square) with the 
cavity measurement on a macroscopic trilayer (open square). 
Dotted curve is a parameter free fit of the FMR using the 
Kittel formula [Eq. ©]. 



as such, since the frequency uo s = 2eVo/h agrees, without 
fitting parameters, with the Kittel formula [3] 



(7) 



for the in-plane magnetic field dependence of the uni- 
form FMR mode when the anisotropy field is perpendic- 
ular to the plane. The anisotropy field Hk = 4900 G 
and the magnetisation at saturation Ms = 930 G are 
both determined directly from the static magnetisation 
data, and j e = /ie/ft, where \±b is the Bohr magne- 
ton. For comparison, the ferromagnetic resonance of a 
macroscopic Nb/PdNi/Nb trilayer has been measured in 
a conventional 9.5 GHz cavity spectrometer at 10 K with 
field applied parallel to the substrate. The cavity FMR, 
shown in the bottom insert of Fig. 3, occurs at 2160 G, 
again exactly as predicted by Eq. 0. Displayed in the 
top insert of Fig. 3 is the comparison of the resonant 
mode in the Josephson junction (solid square) and in the 
macroscopic trilayer (open square). The dotted curve 
shows the frequency of the FMR mode calculated from 
the Kittel formula [Eq. (j7|)]. The spectra presented in 
the main part of Fig. 3 contain an extrinsic broaden- 
ing caused by a lock-in modulation voltage of ~ 1 /iV. 
For the ac modulation voltage of 0.5 /iV, the junction 
resonance width saturates at 0.5 /iV, which corresponds 
to the conventional resonance width (150 G). The sig- 
nal amplitude corresponds to a resonant susceptibility of 



approximately 10, consistent with the FMR mode mea- 
sured in a microwave cavity and reported in the bottom 
insert of Fig. 3. 
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FIG. 4: (color online), (a) The dynamical resistance of the 
ferromagnetic Josephson junction with applied external mi- 
crowaves. Pronounced dip is the Shapiro resonance, and ar- 
rows indicate sideband resonances at the same frequency as 
the n — 2 mode from Fig. 3. The external microwave frequency 
is 17.35 GHz and the temperature 35 mK. (The solid curves 
are the measured derivative of the dynamical resistance from 
Fig. 3, and dotted curves are a fit [Eq. JB])]. The spin wave 
frequency increases with field. Insert: The field dispersion 
relation of the modes. The solid curve is Eq. (7) when the 
spatial dependence of the FMR modes is taken into account. 

In order to demonstrate that the magnetic system is 
coupled to the super- but not to the normal current, 
we have performed Shapiro step [2] measurements, re- 
ported as the dynamical resistance dV/dl in Fig. 4(a). 
The junction is irradiated with microwaves of frequency 
v = 17.35 GHz at 35 mK. The Shapiro steps arise from 
the mixing of the microwave signal with the ac Josephson 
effect and are smaller replicas of the zero- voltage current 
step displaced from zero voltage by V n (h/2e)(j, where 
cj = 27t^, and n is an integer. We do not observe half in- 
teger Shapiro steps, indicating negligible higher harmon- 
ics in the current-phase relation. However, as expected 
within the RSJ model, the ferromagnetic resonance can 
be exited at voltage V ns = (h/2e)uo s /n The sub- 
harmonic for n = 2 is visible in the spectrum in Fig. 3. 
As shown in Fig. 4(a), it is reproduced as a side-band to 
each regular step when V = (Ti/2e)(nuo ± uo s /2). Exper- 
imentally, we do not have available a high enough fre- 
quency to separate similar side-bands for the main FMR 
mode at uo s . 

Finally, the field dependence of the resonance at uo s /2 
has been studied in more detail in the second deriva- 
tive, d 2 V/dI 2 [Fig. 4(b)], where the minima correspond 
to Vis — (2e/h) uo s /2. Measurements were limited in field 
due to the rapid decrease of the critical current above 
800 G. 

In the insert of Fig. 4(b), we show (1 / '2tt) u s / '2 as a 
function of the applied magnetic field. The error bars 



are due to the drift of the amplifier. The solid curve 
is Eq. (|7j), without fitting parameters, with the spatial 
dependence of FMR taken into account. The spatial de- 
pendence of the spin-waves leads to an additional term 
to Eq. Q given by a/c 2 , where k = (7rd/&o)H is the spin 
wave momentum and a = E ex b 2 , where E ex = 50 meV 
is the PdNi exchange energy and 6 = 0.1 nm the lattice 
constant. Since the width of the junction is only about 
500 nm, this leads to a small but finite correction to the 
uniform FMR energy which is larger than the direct ef- 
fect of the applied dc field. Illustrated in this manner 
is the direct determination of spin- wave dispersion using 
the present technique. 

In conclusion, we have demonstrated the dynamical 
coupling of the superconducting phase with the spin 
waves in a ferromagnet and measured their dispersion. 
We have performed a photon free FMR experiment on 
about 10 7 Ni atoms, which would be infeasible with 
standard FMR techniques, and have illustrated a new 
methodology for the study of spin dynamics. There are 
direct and implied applications to spintronics and nano- 
magnetism [15j . 
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